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I.  Introduction 


Background 

The  concept  and  the  requirements  for  safety,  durability  and  man- 
agement cf  military  aircraft  structures  are  summarized  by  Coffin  and 
Tiffany  (Ref  1) . As  reflected  there,  the  durability  analysis  as  well 
as  scheduled  maintenance  for  the  structure  are  based  on  some  assumed 
usage  of  the  fleet.  Usually  the  actual  usage  of  an  individual  air- 
craft and  the  fleet  will  differ  from  the  usage  for  which  the  aircraft 
was  designed.  In  order  to  provide  updated  usage  data  for  force  manage- 
ment, MIL  STD  1530A(11)  (Ref  2:23-25)  requires  establishing  aircraft 
structures  tracking  programs . Such  a program  provides  baseline  oper- 
ational load  spectra  for  average  fleet  usage  as  well  as  an  Individual 
Airplane  Tracking  (IAT)  program,  which  provides  crack  growth  information 
concerning  each  individual  aircraft.  The  IAT  is  the  feedback  link 
necessary  for  the  scheduling  of  structural  maintenance. 

A typical,  current  technology,  IAT  program  is  that  of  the  F-4 
aircraft.  This  program  was  developed  as  part  of  the  Damage  Tolerance 
Assessment  Study  for  F/RF-4  C/D  and  F-4E(S)  aircraft  (Ref  3,  4).  The 
program  uses  as  tracking  devices  counting  accelerometers  installed  on 
each  aircraft  and  flight  data  recorders  installed  on  12%  of  the  fleet. 

A crack  growth  analysis  is  performed  for  critical  structural  locations. 
Stress  exceedances  spectra  for  each  critical  location  and  for  different 
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usage  categories  were  developed  from  the  flight  recorded  data.  These 
spectra  along  with  counting  accelerometer  exceedances  readings  are 
used  to  determine  the  stress  history  at  critical  locations.  Using 
these,  a crack  growth  analysis  for  an  individual  aircraft  may  be  per- 
formed . 

This  tracking  system,  described  in  more  detail  in  other  places 
(Ref  3,  4,  5,  6:64-66),  is  very  indirect  in  defining  stresses  and  in- 
volves many  sources  for  error.  Further,  the  enormous  amount  of  data 
to  be  collected  in  order  to  get  the  stresses  adds  to  program  cost.  A 
device,  the  crack  gage,  has  been  proposed  which  might  provide  a more 
direct  and  efficient  technique  for  use  in  IAT  programs.  It  has  been 
estimated  that  using  crack  gages  in  the  F-4  IAT  program  might  save 
7.7  million  dollars  in  five  years  (Ref  6;79). 

Crack  Gage  Concept 

A preliminary  concept  of  a crack  gage  was  patented  by  Smith  (Ref  7) 
in  1976.  A more  practical  concept  for  using  a crack  gage  for  monitoring 
structure  crack  was  introduced  by  Crane,  Gallagher  and  Grant  (Ref  8). 
Figure  1 shows  the  concept  of  the  crack  gage . It  consists  of  a pre- 
cracked thin  metallic  sheet  with  ends  bonded  to  a structure.  Since  the 
structure  and  the  gage  will  be  exposed  to  the  same  load  history,  it  is 
possible  to  relate  the  crack  growth  in  the  structure  to  that  in  the 
gage.  The  readily  observed  crack  growth  in  the  gage  can  then  be  used 
as  a prediction  of  crar’k  growth  in  the  structure . 

Efforts  in  various  directions  were  made  in  evaluating  the  crack 
gage  concept  since  it  has  been  introduced.  Torvik  (Ref  9,  10)  developed 
an  analytical  solution  for  the  stress  intensity  factor  for  displacement 
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Center  Cracked 
Uniform  Thickness 


controlled,  edge  cracked  gages,  Ori  (Ref  11)  and  Ashbaugh  and  Grandt 
(Ref  12)  evaluated  the  concept  using  both  edge  cracked  and  center 
cracked  gages . 

Stepped  center  cracked  gages  are  now  being  investigated  under 
two  U.S.  Air  Force  contracts,  one  with  Boeing  Wichita  (Ref  13)  and  one 
with  McDonnell  Douglas  Corporation  (Ref  14) . The  contract  with  McDonnell 
Douglas  will  provide  evaluation  of  a center  cracked  stepped  gage  on  a 
F-4  aircraft  during  a full  scale  fatigue  test. 

Statement  and  Scope  of  Problem 

It  has  been  observed  (Ref  9,  10)  that  relatively  high  stress  inten- 
sity factors  for  an  edge  cracked  rectangular  gage  may  be  obtained  if 
the  gage  is  long,  but  the  magnitude  then  varies  considerably  as  the 
crack  grows.  On  the  other  hand,  constant  stress  intensity  factor  vs 
crack  length  may  be  obtained  for  relatively  short  gages,  but  the  gage 
sensitivity  is  then  very  low.  For  analytical  purposes,  such  as  pre- 
dicting the  remaining  life  of  a structure,  a constant  stress  intensity 
factor  vs  crack  length  is  desired.  It  was  suggested  that  an  edge-cracked 
gage  of  a trapezoidal  planform  might  satisfy  the  requirement  for  a 
higher,  but  constant  Kj  vs.  crack  length,  since  the  "effective"  length 
of  the  gage  could  vary  as  the  crack  grows.  The  stepped  gage  was  also 
considered  since  the  thicker  section  serves  as  an  amplifier  of  the 
strain  in  the  cracked  section  from  the  far  field  value . The  gain  may 
be  controlled  by  varying  the  thickness  ratios  and  length  ratios  of 
the  segments  of  the  gage . 

The  gages  that  had  been  developed  for  the  activities  previously 
mentioned  were  basically  developed  empirically.  It  became  apparent 
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that  an  analytical  method  that  could  be  used  to  obtain  stress  intensity 
factor  solution  for  stepped  center  cracked  gages  was  required.  Gages 
with  different  stress  intensity  factor  sensitivities  are  required  so 
that  the  gage  can  be  matched  to  the  specific  usage , i.e.,  applying 
gages  of  very  high  sensitivity  to  structures  to  be  subjected  to  a large 
number  of  cycles . 

In  this  study , a means  of  obtaining  analytical  predictions  of  the 
stress  intensity  factor  for  the  two  geometries , trapezoidal  planform 
and  stepped  gages , is  developed.  These  methods  are  applied  to  the  design 
of  crack  gage  and  the  results  for  the  stepped  gage  were  experimentally 
verified. 
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II . Stress  Intensity  Factor  Solution  for  Cracked 
Flat  Sheet  Subjected  to  Mixed  Boundary  Conditions 

Introduction 

A method  utilizing  Williams1  (Ref  15,  16)  stress  functions  along 
with  energy  principles  was  developed  by  Torvik  (Ref  9,  10)  to  obtain 
a stress  intensity  factor  solution  for  a finite,  cracked  elastic  sheet. 
Solutions  were  provided  for  a rectangular  edge  cracked  sheet  subjected 
to  mixed  boundary  conditions.  Part  of  the  boundary  was  subjected  to 
prescribed  tractions  as  shown  in  Figure  2A , herein  referred  to  as  the 
type'  1-mixed  boundary  condition.  The  solution  for  a center  cracked 
sheet,  which  is  the  center  section  of  the  stepped  cracked  gage,  led  to 
a different  mixed  boundary  condition,  herein  referred  to  as  type  2, 
and  shown  in  Figure  2B . For  this  problem,  a part  of  the  region  of  the 
physical  surface  is  subjected  simultaneously  to  prescribed  tractions 
as  well  as  to  displacements.  A brief  outline  of  the  method  developed 
by  Torvik  (Ref  9,  10) , for  the  type  1-mixed  boundary  condition  problem, 
and  the  necessary  extension  for  the  type  2-boundary  condition  problem, 
will  be  presented  here. 

Method  of  Solution 

Consider  a thin  elastic  sneet  for  which  plane  stress  conditions 
hold,  with  a crack,  starting  at  an  edge.  A coordinate  system  is  set 
with  origin  placed  at  crack  tip.  The  crack  faces  are  located  at  0=±tt 
0£r£a,  and  are  free  of  tractions.  Williams  (Ref  15,  16)  has  shown 
that  for  an  infinite  domain,  stress  fields  can  be  obtained  by  using  a 
stress  function  of  the  form 

X(r,  e,  X)  = rA+1  F(6,  X)  (1) 
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A = m/2  for  m = 1,  2,  3,  ... 


(2) 


The  stress  function  may  be  divided  into  even  and  odd  parts.  Since 
the  geometries  of  interest  for  this  study  are  symmetric  in  respect  to 
X axis,  only  the  even  part  of  the  stress  function  is  of  interest.  In 
addition,  displacements  corresponding  to  the  stress  function  may  be 
found  from 

2pur  - - + <1  - s)r||  (3) 


2yU0  = - - lX  + (1  - s)r2 


r 00 


3r , 


(4) 


y = shear  modulus 


s = v/l+V 


V = Poisson  ratio 


2 3 . dib. 

V2x  = -“-(r— *-) 
3r  30 


(5) 


The  even  parts  of  X and  ^ * as  well  as  the  corresponding 
stresses  and  displacements,  are  shown  in  Appendix  A.  Each  stress  and 
displacement  is  made  up  from  an  even  and  an  odd  part.  As  stated  by 
Torvik  (Ref  9:3-4,  10),  the  stresses  a and  the  displacements  U may 
be  represented  in  the  form: 
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where  each  am  satisfies  field  equations  and  the  traction  free  require- 
ment on  the  crack  face , and 
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Each  of  Um  satisfies  strain  displacement  equations  and  stress 
strain  laws,  and  the  arbitrary  coefficients  d^  are  related  to  the 
stress  and  displacement  functions  coefficients  presented  in  Appendix 
A as  Am  through  the  relation 

di  = A^ , = -A2  f d^  = -A^ , d^  = A4 , dg  = A5 , dg  = -A^ , etc . 

In  order  to  find  a stress  intensity  factor,  defined  as, 

Kj  = a00(0  = 0 ,r)  } (8) 

r-K) 

it  will  be  required  to  find  the  stress  Gqq  for  a certain  boundary 
condition  problem.  More  specifically,  the  coefficients  in  Eq  (6) 

are  to  be  obtained. 

After  applying  the  limit  procedure,  the  only  term  to  remain  is 
the  one  multiplied  by  d^  . All  other  terms  vanish  since  they  include 
r^  expressions , pCll/2 . Hence , 

Kx  = /2tP  A±  = -v^F  d (9) 

The  solution  of  the  coefficients  for  a mixed  boundary  condition  of 
type  1 was  given  by  Torvik  (Ref  9:20-21,  10)  utilizing  Reissner's  prin- 
ciple. The  interpretation  of  Reissner's  principle  by  Fung  (Ref  17:299- 
300)  is  that  the  surface  S enclosing  the  domain  V , can  be  divided 
into  Srp  and  Su  , where  tractions  are  prescribed  on  ST  and  displace- 
ments on  Su  . A more  general  interpretation  will  be  presented  here. 

In  the  type  2 mixed  boundary  condition,  a surface  region  on 

the  surface  S can  be  subjected  simultaneously  to  prescribed  tractions 
and  displacements  as  well.  Hence,  it  may  be  considered  a subset  of  ST 
as  well  as  Su  , with  the  relations : 
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s = ST  + su 


ST  - CT  + CTu 


Su  cu  + ^Tu 


(10) 

(11) 

(12) 


where  Gp  is  the  region  where  only  tractions  are  prescribed,  and  Cu 
is  the  region  where  only  displacements  are  prescribed. 

Consider  the  functional  (essentially  Reissner's  function) 

J " Jt  W(e±j)  - F-Ujdv  - 1/2  (Ui;j  + Ujf.)]dv  - 


- I Ti*  - J a±j  v.  (u±  - u±*) 

* r>  ' r> 


dS 


(13) 


u 


where  stresses , displacements  and  strains  are  to  be  independent  vari- 
ables, and  0^=0^.  The  elastic  equilibrium  will  be  obtained  by  setting 
J - 0,  or  J will  be  stationary  with  J = 0 „ The  are  introduced  as 
Lagrange  multipliers . 

Sj  = 0 = 6e.  . - Fj  6u • ) dv  - So  . . [e.  . - 1/2(U.  . + U.  .)]dv  - 

J 3ef  j ID  1 1 J ID  ID  irD  D#i 

v v 

~f  °if6eij  - V2(6uifj  + 6u.,.)]dv  -J  T^Su.ds  - 

v '"St 


~f  - Ui*)ds  - J v.a.^  u.ds 

* ~ ^ r* 


(14) 


3u 


Rearranging  and  applying  the  divergence  theorem  to  third  integral,  we  find 


6J  = 0 =f  {f-  - aij)6a.jdv  - f (PA  + ^Su.dv  - fso  i;I^ 

V V ^ J T- 

- l/2(Uifj  + U..  jjdv  + dS  -/  Ti 

^ C 


' ~ ei  j - 


- j (Ui  - Ui*)6TjLdS  - j Ti6Ui 


s T 

dS 


(15) 


3U 
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Realizing  that  the  first  three  integrals  represent  the  stress-strain 
law,  equation  of  equilibrium  and  strain  displacement  equation  respec- 
tively, they  may  be  omitted  from  Eq  (15)  if  we  choose  the  stresses  to 
be  represented  by  Eq  (6) , and  displacements  by  Eq  (7) , along  with  the 
corresponding  strain. 

Eq  (15)  reduces  to: 


6u±ds 


♦SU-jdS 


- T.  *)  <5U. dS 
x x 


(U±  - Ui*)6TidS 
u 

- U.*)6T.dS  = 0 
x x 


6U±dS=0 


(16) 


(17) 


Applying  Eqs  (11)  and  (12) 


f (Ti  - Ti*)  6th  dS  -j  (Ui  - U^)  *6T^dS  + 


u 


[(T(±)-  T(i*).<SU,,,  - (U,J4  - U , _.  * ) *6T/jx  |dS  = 0 (18) 


(i)  ' (j)  (j)'  ~ (j) 


]■ 


Tu 


where  j^i,  and  the  parentheses  indicate  suspension  of  the  summation 
convention.  On  , one  component  each  of  displacement  and  traction 

are  given. 

Eq  (18)  may  be  used  to  find  a Kj  solution  for  problems  with  mixed 
boundary  condition  of  type  2 . 

Use  now  the  stresses  and  displacements  as  defined  in  Eq  (6)  and  (7) 
to  find  variations  of  stress  and  displacement. 


M 


d*cr.  . 

P ID 


p=l 


(19) 
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(20) 


Substituting  Eqs  (6) , (7) , (19)  and  (20)  into  Eq  (18) : 
M M 


f tWj  £ - v>  iLsaPuiPlas 

' CT  Ipl  p=l 

' Cu  m=l  P=1 

M M 

f V(«)-  «*"(?)- 

SC  J m=l  ^ P=1 

Tu 


M 

_ i_  y 

2V-  i. 


m 


2y  'fei  2yU(3y  P= 


M 

*>)viZ_.  6dPa-  / 
J>  *=1 


dS  = 0 


(21) 


Since  6 dp  is  arbitrary,  make  the  choice  6dp  = 0 , p?*q,  6 dp  - 1 


if  p=q  and  then 
M 


M r 

(VjZ  Vi/V  &S  - 1 ^n,uim,vjaijqas  - 

SfC*'1 


~ ^m=J 
C t-C 
ur  Tu 


= I T.*»U±qdS  +1  tTi**UiqdS  - 2yUi*Tiq]dS  - 2 J U^-Tj^dS 

Jr  1 ^0  X 1 JC 

T M U (22) 


Eq  (22)  may  be  represented  as : 


where  A^ 


M 


m=l 


F 

q 


is  a symmetric  matrix. 


(23) 
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- / r / . ni  m m cr  i 

V - I I'Vrr  + Vre  ,Ur  + <Ver  + Vee,uelds  ■ 


-I 


c c 

T+  Tu 


+ ve°eeq>  + '’.“‘Vrr'1  + 'V0ldS 


CufCTu 


r r rr  6 r0 


(24) 


The  reader  is  reminded  that  only  one  product  of  force  and  displacement 
is  to  be  evaluated  in  the  contribution  of  to  each  of  these  integrals . 

F<I  = J (T*rUrq  + T0*Ueq)dS  - 2]f  (Ue*TQq  + Ur*Trq)  dS 

JC 


■^Tu 


[Tr*Urq  + T0*U0q 


2y(U*T0q  + U*T^q)]dS 


r r 


(25) 


Since  the  displacements  are  prescribed  elsewhere,  the  crack  tip 
has  to  be  allowed  to  translate.  Thus,  as  stated  (Ref  9:20,  10),  we  may 
introduce  a translation  term  into  the  expansion,  which  will  be 


0 


2viUr 

" Vi  “ 

c1qCos0 

(26) 

to 

B 

CD 

= d0ui  - 

-dQsinG 

(27) 

with  arr°  , aee°  , 

ar0° 

zero. 

By  solving  Eq  (23)  along  with  Eqs  (24)  and  (25)  for  the  coefficients, 
dJn  a solution  for  a finite,  two  dimensional,  cracked  elastic  sheet  with 
mixed  boundary  conditions  of  type  2 is  obtained.  Note,  that  by  elimi- 
nating the  integral  on  C^u  region,  we  get  a type  1 mixed  boundary 
problem,  and  A^  and  Fq  will  have  the  same  form  as  shown  (Ref  9:23,  10) 
previously. 

Since  the  particular  gage  geometries  for  this  study  are  the  trap- 
ezoidal gage  and  the  stepped  gage , the  formulation  for  getting 
and  Fq  for  these  geometries  will  be  shown. 
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The  basic  boundary  conditions  for  the  trapezoidal  gage  are  essen- 


tially as  shown  in  Figure  2a,  anu  the  formulation  was  presented  (Ref  9: 
23,  10) . 

Fq  = 2V/UVCOS6  + VsinO)(Vrarr^VeaerS  + 

+ (-Ux*sin6  + UyXcos8)  • (vrar6q  + v^^ldS  (28) 


When  evaluating  A^  for  the  discussed  geometry,  crack  faces  are  free 
of  traction.  The  first  integral  of  Eq  (24)  will  be  evaluated  on  edged 
1 and  3,  while  second  integral  will  be  evaluated  on  edge  2. 

The  case  with  boundary  conditions  shown  in  Figure  2b  represents  the 
center  portion  of  the  stepped  gage.  Edges  1 and  2 are  part  cf  CT  , 
where  tractions  are  prescribed,  while  edge  3 is  part  cf  Cnpu  , 

Taking  into  account  boundary  condit ions , only  a part  of  the  integral 
CTu  exists  on  edge  3.  Due  to  boundary  condition  presentation,  A^q 
and  will  be  represented  in  cartesian  coordinates,  and  Eq  (24) 

becomes 


- umTq]dS 
x x 


(29) 


With  Cu  = 0 , and  with  the  prescribed  variables  on  CTu  substituted 
into  Eq  (25) , we  get 


■/. 


(T  *U 
x * 


+ T 


*uq) ds 


y y 


(30) 


Those  last  two  equations  will  be  used  to  solve  the  stepped  crack  gage 
stress  intensity  problem. 
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Ill „ Stress  Intensity  Factor  Solution  for  A Trapezoidal , 

Edge  Cracked  Crack  Gage  with  Prescribed  End  Displacements 

The  edge  cracked  rectangular  gage  has  a tendency  to  give  low  stress 
intensity  factor  values  for  short  gages , while  the  sensitivity  improves 
for  relatively  long  gages  (Ref  9:34-46;  10).  This  has  also  been  ob- 
served in  experiments  (Ref  6:73).  The  disadvantage  of  long  gages  is 
that  Kj  does  not  remain  constant  and  even  decreases  as  crack  grows . 

This  property  is  not  desired  for  a crack  gage.  It  was  thought  that 
giving  the  gage  a trapezoidal  shape  as  shown  in  Figure  3,  a better  con- 
trol on  the  Kj  variation  with  crack  growth,  and  a higher  Kj  might 
be  achieved. 

The  suggested  configuration  consists  of  a thin  gage  bonded  to  a 
structure  along  the  line  AB.  The  structure  is  loaded  by  stress  Os 
and  displacements  are  transferred  to  the  gage  through  line  AB.  The 
stress  intensity  factor  solution  for  this  case  was  obtained  by  using 
the  method  shown  in  Section  II.  Using  Eq  (28)  along  with  Eq  (24) , 

Eq  (23)  gives  a solution  for  the  coefficient  . 

The  prescribed  displacements  are 

* ^s 

u*  - — y (3D 

* Es 

U*  = -XV  — (32) 

Es 

where  y = Ha  - tgcr  (b  - x)  and  v=.3  . (33) 

The  matrix  was  evaluated  by  a computer  program  (see  Appendix 
B)  using  an  expansion  of  thirty-six  coefficients,  including  those  for 
rigid  body  displacements . The  integration  was  performed  over  300 
boundary  points.  The  variables  for  this  study  were  chosen  to  be  the 
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angle  a , which  was  assumed  to  control  Kj  changes  as  the  crack  grows , 
and  also  to  control  the  magnitude  of  Kj  , since  the  gage  length  dimen- 
sions are  changing.  Former  results  for  a rectangular  gage  (Ref  9:34, 

10)  indicated  that  a gage  with  an  aspect  ratio  of  2 Ha/b=2  might  be  a 
good  one  for  the  study,  since  for  that  gage  was  fairly  constant 

as  crack  length,  a , increased. 

Results  and  Discussion 

Non-dimensional  stress  intensity  factor  results  for  Hc*/b=l  are 
presented  in  Figure  4 . It  can  be  observed  that  changes  in  a in  the 
range  of  160°  - 200°  have  a very  little  influence  on  the  magnitude  of 
Kj  . The  general  behavior  of  Kj.  vs . a improves  a little  and  for 
a = 200°,  Kj/asv^  is  almost  constant  on  the  range  .4<a/b<.8. 

Results  for  Ha/b  - 1 are  presented  in  Figure  5.  It  can  be  observed 
that  in  comparison  to  former  results  the  changes  in  Kj  are  relatively 
large  as  a changes  from  165°  to  195°.  This  can  be  related  to  the 
fact  that  under  the  displacement  loading,  the  will  be  more  sensitive 
to  length  changes  on  Hc  than  on  Ha  since  the  former  directly  affects 
the  crack  opening,  and  Kj  may  be  considered  proportional  to  crack 
opening . 
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Figure  5.  Stress  Intensity  Factor  for  Trapezoid  Crack 

Gage  Ha/b  = 1. 
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IV.  Stress  Intensity  Factor  Solution  for  a_  Stepped 

Center  Crack  - Crack  Gage , with  Prescribed  End  Displacements 

The  stepped  center  crack  gage  concept  is  described  in  Figure  6. 

The  gage  consists  basically  of  three  segments.  The  center  segment  of 
reduced  thickness  which  includes  the  crack,  and  two  thicker  end  seg- 
ments . 

Part  of  each  thicker  segment  is  bonded  to  the  structure,  and  a 
part  (of  length  L0)  is  free.  It  is  assumed  that  the  gage  is  clamped 
to  the  structure  along  lines  AB  and  GH.  This  assumption  is  only  an 
approximation  since  the  actual  load  transfer  from  the  structure  to  the 
gage  should  be  modeled  as  loads  transferred  through  springs,  where  the 
adhesive  layer  is  the  spring  element.  However,  the  exact  model  wTas 
beyond  the  scope  of  this  study.  Displacements  of  line  AB  relative  to 
line  0-0  are  assumed  to  be  transferred  completely  to  the  gage  with  no 
losses  in  the  bond  line,  i.a.r  the  bond  line  does  not  rotate. 

Due  to  the  symmetry,  the  center  segment  of  the  gage  needs  to  be 
solved  only  for  a quarter  of  the  plate,  as  shown  in  Figure  7.  The 
entire  outer  segment,  shown  in  Figure  8,  was  analyzed  in  order  to  obtain 
better  accuracy  with  the  finite  element  analysis. 

Methods  of  Solution 

Two  analytical  methods  will  be  used  to  obtain  the  complete  solution. 
The  Finite  Element  Method  will  be  applied  to  solve  forces  and  displace- 
ments of  the  outer  segment  of  the  gage  and  the  method  described  in 
Chapter  II  will  be  used  to  solve  for  stress  intensity  factors,  stresses 
and  displacements  for  the  center  segment.  Bowie  et  al  (Ref  18:767-772) 
introduced  also  the  idea  of  getting  stress  intensity  factors  by  using  a 
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Figure  6.  Center  Cracked,  Stepped  Gage  Geometry 
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stitching  technique.  They  used  one  analytical  method  to  solve  stresses 
and  displacements  for  cracked  and  uncracked  part.  By  stitching  the 
solutions  stress  intensity  factors  were  found  for  center  cracked  uni- 
form thickness  plate  subjected  to  uniform  stress  loading. 

Since  the  line  CD  is  actually  the  connecting  line  between  the 
two  segments,  we  will  demand  displacements  and  forces  produced  by  the 
two  methods  to  be  equal  on  this  line.  Let  the  line  CD  be  discretized 
into  2n  elements,  each  of  length  Dx,  Each  element  will  be  represented 
by  a nodal  point,  considered  to  be  at  the  center  of  the  segment,  at 
which  forces  and  displacements  are  assumed  to  be  applied. 

The  outline  for  the  solution  is  as  follows: 

{u*}  = {uc}  + {uQ}  (34) 

where  U*  , Uc  are  the  prescribed  relative  displacement  of  line  AB 
(Figure  6)  and  relative  displacement  of  line  CD  with  respect  to  inertial 
frame  x,  y,  respectively.  Uq  is  the  relative  displacement  of  line  AB 
with  respect  to  CD.  Using  flexibility  matrix  of  center  segment.  [sc]  , 
and  flexibility  matrix  of  outer  segment  [sQ]  , we  have 


{uc}  = 

tsc] 

{pc} 

(35) 

{U  } - 

oJ 

[s0J 

Oo> 

(36) 

On  line  CD,  the  forces 

{p0}  : 

■ tPc) 

(37) 

After  substitution 

{u*}  - [fs0]  4 

. rs  ] 1 
LScJj 

i (Pc) 

(38) 

Aid  the  forces  and  displacements  are  expressed  by  x and  y direction 


components  consequent Iv  all  matrices  ere  partitioned  and  the  number  of 


the 


of  nodal  points  under  consideration . 


rows  or  columns  is 


number 


Suppose  a unit  load  is  applied  to  each  element  of  the  center  seg- 
ment, successively.  Using  the  technique  described  in  Chapter  II,  the 
coefficients  of  the  expansion  may  be  solved  for  each  loading  case. 

From  the  coefficients  solved  for  a unit  load  at  nodal  point  k,  the 
first  is  of  particular  interest  and  will  give  the  stress  intensity 
factor  for  unit  load  at  nodal  point  k . When  we  repeat  the  procedure 
for  all  nodal  points,  we  may  construct  an  array  {d}  whose  k'th  element 
is  the  contribution  of  load  at  point  k to  the  stress  intensity  factor. 
The  stress  intensity  factor  for  the  center  segment  under  uniform  dis- 
cretized unit  load  is  obtained  by  summing  the  elements  of  {d}.  If  the 
load  at  each  nodal  point  is  different,  the  stress  intensity  factor 
solution  may  be  written  as : 


Kj  = 2tt{d}T  {Pc  } 


(39) 


or, 


Kj  = 2tt{d}T  [ [Sc  ]+  [SQ  ] ]"1{U*} 


(40) 


Thus,  the  solution  for  the  stepped  gage  with  prescribed  end  dis- 
placements requires  that  the  flexibilities  [ Sc  ] and  [sQ]  be  determined, 
A prescribed  stress  T*  , on  k'th  element,  Dx  long,  may  be  con- 
verted into  a "concentrated"  stress  at  the  nodal  point  xk  producing 
a unit  force  per  unit  thickness  by  defining  a loading  function. 

T*  = 6(x  - Xfc)  T*  = 6(x  - x^)  (41) 


so  as  to  produce  a unit  force  in  each  coordinate  direction  on  any 
element,  centered  on  xk 
xk+Dx/2 


6 (x  - xk)dx  - 1 


>:y^Dx/2 


(42) 
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For  simplicity  a unit  thickness  is  considered.  The  vector 
can  be  computed  by  Eq  (30) , with  displacements  evaluated  at  x^  . 

The  matrix  A^  can  be  evaluated  by  Eq  (24)  and  the  system  of  equations 
can  be  solved  for  the  coefficients  d^  by  using  Eq  (23) . With  the 
set  of  coefficients  known,  the  actual  displacements  at  every  other 
nodal  point  1 can  be  evaluated.  By  repeating  the  procedure  for 
k = l...n,  we  contruct  the  flexibility  matrix  [3C] , along  with  the 
array  {d}  which  is  constructed  from  the  coefficients  d^  for  each  of 
k load  points.  The  flexibility  matrix  [Sc]  is  actually  factored  such 
as  to  account  for  actual  thickness  ratio  between  center  and  outer  seg- 
ments of  gage.  The  outer  segment  flexiblity  matrix  is  contructed  by 
using  "ANALYZE"  (Ref  19) , a finite  element  analysis  computer  program 
developed  at  the  Air  Force  Flight  Dynamics  Laboratory.  The  matrix  was 
established  by  constraining  one  end  of  the  segment  completely  to  simu- 
late the  bond  line,  and  the  other  end  was  loaded  with  unit  forces  at 
discrete  points . 

Numerical  Evaluation 

The  problem  introduced  for  the  center  segment  was  solved  ty  com- 
puter programs  which  are  described  in  Appendix  B.  The  objective  of  the 
program  was  to  evaluate  the  matrix  A^  , given  in  Eq  (29) , to  solve 
Eq  (23)  for  the  coefficients,  and  to  construct  the  flexibility  matrix 
[sc]  in  the  manner  previously  described. 

The  matrix  A^  was  evaluated  by  a numerical  integration  along 
edges  1,  2,  3,  using  Simpsons'  third  rule.  The  boundary  was  divided 
into  300-400  points.  A slight  sensitivity  to  the  number  of  integration 
points  was  observed  (for  short,  Lc)  when  elements  of  different  lengths 
were  used  to  perform  the  integration  on  separate  edges . 
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The  number  of  terms,  M,  used  for  the  expansion  was  18.  The  number 
of  nodal  points  used  on  edge  2 was  20.  Consequently,  the  flexibility 
matrixes  were  dimensioned  40x40,  partitioned  into  x and  y components. 

The  flexibility  matrix  of  the  outer  segment  was  developed  for  a 
whole  outer  plate,  and  the  corresponding  40x40  matrix,  £s0J  , was 
extracted. 

The  prescribed  displacement  U*  was  chosen  to  be  1 in  longitudinal 
y direction  and  zero  in  x direction.  As  has  been  shown,  releasing 
the  transverse  constraints  would  have  little  effect  on  results  (Ref  9: 

32,  10),  since  the  stress  intensity  factor  is  driven  primarily  by 
displacments  perpendicular  to  crack  plane. 

Stress  intensity  factors  were  computed  for  center  segments  having 
lengths : 

Lc  = .4;  .6;  .8;  1.0;  1.25 

The  outer  segments  having  lengths 

Lq  53  1.0;  1.267;  1.6;  2.0;  4.0 

and  for  the  thickness  ratios 

tc/tc  = .2;  .5;  .75 

The  width  in  all  cases  was  taken  to  be  b=2.0. 

Stress  intensity  factors  were  not  calculated  for  all  combinations  of 
variables . 

Results  and  Discussion 

A sample  of  the  results  of  a convergence  study  for  the  stress  in- 
tensity factor  solution  of  the  center  segment  with  12,  18  and  24  terms 
is  presented  in  Table  I.  No  significant  changes  for  more  than  18 
terms  were  observed.  The  number  of  nodal  points  N,  was  varied  through 
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10,  20 , 30  , 40  and  results  for  stress  intensity  factors  are  compared 
in  Table  II  for  a/b  = .2;  ,6.  A rapid  convergence  may  be  observed 
for  n >10.  Using  results  for  40  nodes  as  a reference,  10  nodal  points 
give  about  63%,  and  20  nodal  points  give  94%  of  the  reference  value 
for  a/b  = .2.  For  a/b  = .6  the  results  for  n = 20  is  98%  of  the 
result  obtained  using  n - 30*  The  convergence  for  longer  cracks  seems 
to  be  better,  since  the  relative  error  of  load  fraction  distributed 
over  the  crack  (governed  by  element  size)  is  smaller  for  long  cracks. 

The  stress  intensity  factor  solution  obtained  for  the  uniform 
discretized  stress  loading  on  the  central  segment  is  given  in  Table 
III,  and  is  compared  to  published  results  for  a uniform  stress  loading 
(Ref  20:11).  Table  III  shows  that  present  results  are  low  by  five 
percent , a factor  that  may  be  attributed  to  the  number  of  nodal  points . 

In  order  to  check  the  overall  approach  cases  of  tc/tQ  - 1 were 
run.  This  case  is  actually  a case  of  a center  crack  gage  with  uniform 
thickness  and  half  length  of  Lc  + LQ,  subjected  to  prescribed  displace- 
ments at  the  end.  Comparison  to  published  results  (Ref  20:16)  indicated 
that  present  results  are  within  10'15%  higher  than  those  published  for  the 
same  aspect,  ratios,  and  uniform  thickness. 

Due  to  the  fact  that  the  gage  is  constructed  from  relatively  short 
plates,  the  finite  elements  solution  w as  checked  thoroughly,  and  re- 
sults were  found  reasonable  (see  Appendix  C) . However,  the  finite  ele- 
ment analysis  was  performed  on  a complete  outer  segment  instead  of  a 
half  since  the  displacements  results  were  more  reasonable  in  this  con-™ 
figurati  on . 

During  the  center  segment  analysis,  attempts  to  invert  the  flexi- 
bility matrix  were  made  and  difficulties  encountered.  However,  no 
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Table  I 


Convergence  of  Stress  Intensity  Factor 

Ki/os^ 


Lq  — .4 


M**  = 12 


1.2620 


b = 2.0 


M - 18 
1.2660 


n = 20 


M = 24 
1.2662 


*n  = number  of  nodal  points  on  half  width 

**M  = number  of  terms  in  expansion 

2M  = number  of  coefficients  in  expansion 


Table  II 

Convergence  of  Stress  Intensity  Factor 
Kj/CTgtd? 


Lc  - .4 

b = 2.0 

M = 18 

a/b 

n - 10 

n - 20 

n = 30 

n = 4C 

.2 

.8497 

1,266 

1.322 

1.30 

.6 

6.607 

6.731 
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Table  III 


Stress  Intensity  Factor  for  Stress  Loading  - Comparison 

with  Published  Results 


a/b 

Lc/b  = .5 

Kt/CT  Zb" 

I'  s 

Present  Results 
Discrete  load 

Rooke  (Ref  20) 
(Data  transferred  to 
same  nondimensional 
form) 

.2 

.84 

• 

00 

CO 

.4 

1.70 

1.78 

.6 

3.15 

3.28 

explanation  for  this  difficulty  was  found.  The  final  results,  as  well 
as  the  intermediate  results  of  the  analysis,  were  within  reasonable 
agreement  with  other  published  results.  Thus,  the  flexibility  matrix 
was  accepted  in  spite  of  the  unexplained  difficulty  in  computing  an 
inverse . 

Figures  9 through  19  show  the  stress  intensity  factor  results 
obtained  in  this  study.  One  of  the  gages  (Lc/b  = .2  Lc/b  = .63 
tc/t0  = *5)  was  developed  by  McDonnell  Douglas  Corporation  as  part 
of  the  evaluation  of  the  crack  gage  cn  F-4  full  scale  fatigue  test 
(Ref  14) . Another  gage  (Lc/b  = .625  L0 /b  = .625  tc/t0  “ .2)  was 
developed  by  Boeing  Wichita  as  part  of  contract  F33615-77~C-5073  (Ref 
13) . The  results  are  given  in  a nondimensionalized  form  (Kj/a5^E) 
(Es/Eg)  vs.  a/b,  where  as  is  the  stress  of  the  carrying  structure 
at  the  bond  line , such  that  a longitudinal  displacement 

U*  = Lrp  (43) 

E 


30 


is  transferred  to  the  gage,  b is  the  width  of  the  gage  and  a is 
the  length  of  the  crack.  Figure  9 shows  the  typical  benefits  of  the 
stepped  gage  compared  to  a gage  with  uniform  thickness.  For  the  same 
aspect  ratio,  i.e.,  1.5,  the  stepped  gage  gives  a 40%  higher  stress 
intensity  factor  and  is  more  constant.  Furthermore,  the  range  of 
stress  intensity  factor  available  to  choose  for  a specific  need  with 
geometrical  restrictions  is  very  large.  Figures  10  through  18  show 
that  stress  intensity  factor  will  grow  as  thickness  ratio  will  de- 
crease, and  as  L0  grows  for  a given  Lc  . Figure  19  is  typical, 
showing  the  effect  of  change  in  thickness  ratio  change  while  the 
length  ratio  is  held  fixed.  It  may  be  observed  that  changes  in  length 
ratios  of  center  and  outer  segment  appear  to  change  the  magnitude  and 
the  shape  of  SIF  vs.  a/b,  while  thickness  changes  will  affect  primarily 
the  magnitude  only.  It  was  found  that  for  a thickness  ratio  less  than 
tc/t0  = .1,  the  stress  intensity  solution  is  as  if  the  prescribed  dis- 
placement were  applied  directly  to  the  center  segment. 

For  fast  approximations  of  the  stress  intensity  factor,  the  following 
equation  may  be  used  (see  Appendix  D for  details) . 


K / E_ 


I 

ym 

5 \ 


Lo  + LC 


SE 


,£ 

' K b 


£o.  —C 


v'S-a)  + /(l-V^)L 


(44) 


In  this  equation,  Rice's  (Ref  22)  limit  values  for  Kj  in  a 
short  infinite  strip  with  a semi-infinite  crack  is  encountered.  The 
values  of  Kj  obtained  by  this  approximation  are  within  10-15%  lower 
than  the  actual  solution,  for  a/b  > .6.  The  approximation  was  verified 
for  values  of  tc/t0  and  LC/LQ  in  the  range  investigated,  and  is  suggested 
as  a rapid  means  of  obtaining  a trial  design  with  desired  stress  intensity 
factor . 
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a/b  - CP  ACX  LENGTH 

Figure  9.  Comparison  of  Stress  Intensity  Values  for  Center  Cracked 
Uniform  Thi ckness  Gage  and  Stepped  Gage 
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a/b  - CRACK  LENGTH 


Figure  10.  Stress  Intensity  Factor  Results  for  Gage 
with  tc/t0  = . 2 , Lc/b  - . 2 
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# “ V*  " • 

+ _ L0/b  = .63 

© - L0/b  = .8 

© - L^b  = 1,0 
X - LQ/h  = 2.0 


.2  .3  .4  .5  .6  .7  .8 

a/b  - CRACK  LENGTH 

Figure  12.  Stress  Intensity  Factor  Results  for  Gage 
with  t„/tQ  = .2,  Lc/b  = .625 
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.2  .3  .4  .5  .6  .7  .8 

a/b  - CRACK  LENGTH 

Figure  14.  Stress  Intensity  Factor  Results  for  Gages 
with  tc/tG  = .5,  Various  Length  Ratio 
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a) • (s^0/1*) 


1.0  - 


.2 


.3 
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.4  .5 


.6 


.7 


.8 


a/b  - CRACK  LENGTH 


Figure  19 . 


Stress  Intensity  Factor  Results  for 
Lc/Lc  = .315  Lc/Lo  .312 
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V.  Experimental  Results 


The  methodology  for  the  verification  of  analytical  results  of 
stress  intensity  factor  for  center  cracked  stepped  gage  is  based  on 
Paris'  law. 

da/dN  = C(AK)m  (45) 

where  a is  the  crack  length;  N the  number  of  loading  cycles;  C , m 
the  material  constants  and  AK  is  the  difference  of  stress  intensity 
factor  for  each  loading  cycle.  If  we  measure  da/dN  in  an  experiment, 
we  can  find  AK  , provided  the  material  constants  C and  m are 
known.  When  conducting  a constant  amplitude  stress  cyclic  test  with 
R=0  where  R is  amin/amax  then  AK  becomes  actually  the  stress  inten- 
sity factor.  This  sort  of  test  is  to  be  designed  such  that  no  retar- 
dation effects  will  be  experienced  during  the  test.  A series  cf  experi- 
ments were  designed  to  check  the  analytical  results  of  this  study  and 
in  addition  were  a part  of  the  evaluation  of  the  crack  gage  designed 
by  McDonnell  Douglas  for  the  F-4  full  scale  fatigue  test.  Due  to 
technical  difficulties,  this  test  program  is  not  completed  and  inter- 
mediate results  are  not  reported.  However,  Boeing  (Wichita)  as  part 
of  contract  No.  F33615-77-C-5073  (Ref  13)  has  conducted  a series  of 
base  line  tests  and  constant  amplitude  stress  crack  growth  test  for 
a similar  stepped  crack  gage.  Since  their  data  was  applicable  for 
verification  of  the  results  presented  in  this  study,  raw  data  was  ob- 
tained and  reduced  to  compare  to  analytical  results  to  be  presented. 

The  material  used  for  the  Boeing  specimens  was  7075-T651  aluminum. 
The  Paris'  law  constants  we  re  found  to  be  C = ..1479*10  ^ and  m *=  2.49145 
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where  the  units  of  da/dN  and  Ak  are  in/ cycle  and  psi*yin,  respectively. 
Base  line  data  was  fit  to  equation  (45)  with  a computer  program  pro- 
vided by  Dr.  Ashbaugh  and  Dr.  Grandt  from  AFML  (Ref  23) . This  program 
uses  Feddersen's  formula  modified  by  Tada  (Ref  21:2.2)  to  evaluate 
Kj  for  a center  crack  coupon.  The  gage  shown  on  Figure  20,  was  bonded 
to  the  carrying  specimen,  shown  in  Figure  21,  by  FM-73  adhesive.  Con- 
stant load  amplitude  tests  were  conducted  with  Os  = 10,000  psi  and 
R = 0.  Crack  growth  vs.  cycles  was  recorded.  From  a vs.  N data,  da/dN 
was  calculated  and  by  applying  Paris'  law,  Kj  was  determined.  The 
results  shown  in  Figure  22  were  taken  from  two  different  tests  employing 
six  gages , and  where  gages  were  located  back  to  back  on  the  carrying 
specimen.  The  first  test  carrying  specimen,  AFCG-1,  included  two  of 
the  gages.  On  the  second  test  specimen,  AFCG05 , four  gages  were  bonded 
to  the  carrying  specimen.  The  graphical  results  of  Figure  22  show  very 
good  agreement  between  test  and  prediction.  It  has  to  be  noted  that 
there  is  significant  scatter  in  the  experimental  results,  as  is  to 
be  expected.  Scatter  is  inherent  in  this  experimental  procedure  due 
to  variation  in  materials,  loading  accuracy,  crack  length  reading  and 
differentiating  procedure  used  to  obtain  da/dN. 
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(Kj/OgVE)  (Es /Eg)  - Stress  Intensity 


Computed  Line 

AFCG  1 Side  1 Ag2 

AFCG  1 Side  2 Ag2 

AFCG  5 Side  1 Gage  7 

AFCG  5 Side  1 Gage  9 

AFCG  5 Side  2 Gage  7 

AFCG  5 Side  2 Gage  9 


X — o f * 

/ **0^  * 

/«*♦»  ;<* 1 ’ 40  * * 

la  + 4 a 

A +•  A «* 

O * A * 

G° 

■ 

A A 


a/b  - CRACK  LENGTH 

Figure  22.  Comparison  Between  Calculated  Stress  Intensity 
Factor  and  Experimental  Results 
Gage  Dimensions  tc/t0=.2,  Lc/I>o-1.0 
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VI.  Conclusions  and  Recommendations 


This  study  provided  methods  for  obtaining  stress  intensity  factors 
for  trapezoidal  planform  gage  and  stepped  gages.  For  the  stepped  gage 
solution,  a method  was  developed  using  Reissner's  principle  along  with 
finite  element  techniques  which  proved  to  give  satisfactory  results, 
as  verified  by  test.  It  was  clearly  demonstrated  that  the  stepped 
gage  can  be  designed  to  have  stress  intensity  factors  within  a wide 
range . 

Conclusions 

1.  The  trapezoid  shaped  edge  cracked  gage  has  no  significant 
advantage  over  rectangular  gages „ 

2.  The  method  of  solution  introduced  for  stepped  gages  has  been 
shown  to  be  suited  to  the  analysis  of  these  geometries  with  good 
accuracy . 

3.  Tne  stepped  center  crack  gage  provides  a convenient  means  for 
obtaining  desired  stress  intensity  factors  due  to  the  relative  ease 
with  which  substantial  changes  in  K-j-  may  be  obtained  by  varying  the 
geometry,  and  is  to  be  preferred  to  other  geometric  shapes  which  have 
been  considered  for  the  crack  gage. 

Recommendations 

1,  It  is  recommended  that  the  analytical  results  presented  in 
this  study  be  verified  by  further  testing  of  stepped  gages  of  other 
geometries  . 

2.  It  is  recommended  that  the  Kj  solution  method  for  stepped 
gages  introduced  be  included  in  crack  growth  prediction  models  so  as 
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to  provide  an  analytical  model  for  crack  growth  prediction  for 
crack  gages  under  load  sequence  more  complicated  than  constant 
amplitude. 
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Appendix  A 


Stresses  and  Displacements  as  Result  of 
Williams  1 Stress  Function 

For  the  symmetric  type  problems  introduced,  only  the  even  part 
of  Williams'  stress  function  is  of  interest  and  may  be  expressed  as : 

X«  = (-ln""','A2n_irn4l//2{-  cos  (n  - 3/2)0  + — — 3 cos  (n  + 1/2)0}  + 

2n  + 1 


+ (-l)nA„  rn+^{-  cos (n  - 1) 0 + cos (n  + 1)0}  (46) 

2n 


0xr,  0 oq,  U - are  found  by: 


1_  32Xe  i_  9X€ 

r2  W + r 3; 
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are 


r,2 

1 d xe  ^ 1 
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0rQ(l)  " ("1)n”lrn_3/2A2n-l(n  ~ 3/2)*(n  - X/2) * 


{sin(n  + 1/2)0  - sin(n  - 3/2)0} 


(50) 


ar0(2)  = (“l)nrn_1A2n(n) { (n  + l)sin(n  + 1)0  - (n  - l)sin(n  - 1)0}  (51) 


-0  (2) 


ar0  ~ ar0(l)  + ar0(2) 


(52) 
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°rr(l)  = (-1)n"lrn"3/2A2n-l'^(n  “ 3/2)2  “ (n  + 1/2)]*cos(n  - 3/2)6 


(n  - 3/2) (n  - 1/2)  cos (n  + 1/2)6} 


a = (-l)V\*{[(n  - l)2  - (n  + 1) ]cos (n  - 1)6 


(53) 


rr  (2) 


- n(n  + l)cos(n  + 1)6} 


arr  °rr(l)  + arr(2) 
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a00  = (-l)n_1rn"3/2A2  ,(n  + 1/2) (n  - 1/2 ) •{*  + ^ cos (n  + 1/2)0  " 

(1) 


- cos(n  - 3/2)6} 


aOQ,^,  = (-l)nrn-1A0  (n)  (n  + 1)  {cos  (n  + 1)6  - cos(n  - 1)}  (57) 

UU(2)  2n 
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aee  aeed)+  aee(2) 
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Solving  the  equation 

2 3 


v *e  ' 1ST  (r  IT* 


(59) 


get  as  a result 


n-l  rn~3/2  r (n-3/2) 2 - (n  + 1/2) 2 , 
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Substitute  in 
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(62) 


. „ 1 9X  , 2 olp 

2yue  = - 7 96  + (1  ‘ s)r  97 


V = Poisson  ratio  s = V/(l  -i-  V)  , 
y = shear  modules * 

As  a result: 


2pDr(1)  = (-Dn'1rn"1/2A2n_1(  [In  + 1/2)  +£“,•/<”  - 3/2)2  - (n 


- 1/2) ^ 3 


•cos(n  - 3/2)0  - (n  - 3/2)  Geos  (n  + 1/2)0} 
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(64) 


2yur  = 2y(Ur(1)  + Ur(2)) 
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Appendix  B 


Description  of  Computer  Program 

Three  computer  programs  used  to  generate  the  stress  intensity 
factor  in  this  study  will  be  described  here*  Program  COEFDT  was  used 
to  solve  the  trapezoid  case*  COEFINA  was  used  to  provide  the  solution 
for  the  center  segment  of  the  stepped  gage.  Program  TAYLOR  was  used 
to  match  the  segments  of  the  stepped  gage.  A few  subroutines  are 
common  to  both  COEFDT  and  COEFINA.  The  subroutines  and  parts  of  COEFDT 
and  COEFINA  were  developed  by  Torvik  for  a previous  analysis  (Ref  9, 

10)  . 

Subroutine  DPS I 

In  this  subroutine  tangential  displacements  in  Eqs  (66)  , (67)  are 
evaluated.  The  results  are  transferred  to  the  main  program  as  for 
Eq  (66)  and  Eq  (67)  separately. 

Subroutine  DR 

In  this  subroutine  radial  displacements  as  in  Eqs  (63)  , (64)  are 

evaluated. 

Subroutine  SPJR 

This  subroutine  evaluates  radial  stresses  as  in  Eq  (53) , (54) . 

Subroutine  STT 

This  subroutine  evaluates  tangential  stresses  as  in  Eqs  (56) , (57) . 

Subroutine  SRT 

This  subroutine  evaluates  shear  stress  as  in  Eqs  (50) , (51) . 

Program  COEFDT 

This  program  evaluates  the  matrix  A^  as  in  Eq  (24)  . The  eval- 
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uation  is  done  for  each  of  the  three  edges  of  the  three  separately , 
and  the  results  added  to  produce  the  required  matrix  element.  The 
evaluation  is  made  first  on  the  uncracked  edge , and  last  on  cracked 
edge.  Then  Fq  is  evaluated  on  edge  2 via  Eq  (25).  Finally,  the 
system  of  Eq  (23)  is  solved  for  d^,  and  d-j_  is  the  desired  result. 
Program  CQEFINA 

This  program  evaluates  the  matrix  A^  via  Eq  (29)  in  the  same 
order  as  done  in  program  COEFDT.  The  vector  Fq  is  evaluated  on  edge 
2 via  Eq  (30)  by  moving  a unit  load  from  one  nodal  point  to  the  other. 
Fq  is  evaluated  n times  as  the  number  of  nodal  points . Solution 

r , k 

of  Eq  (23)  each  time  gives  finally  an  array  {D}  of  coefficients  d^  . 
In  parallel  to  the  above  described  procedure,  the  program  constructs 
the  flexibility  matrix  [Sc]  for  the  center  segment. 

Program  TAYLOR 

This  program  evaluates  Kj  via  Eq  (40) . 


56 


Finite  Element  Solution 


The  flexibility  matrix  for  the  outer  segment  of  the  gage  was 
obtained  by  using  a slightly  modified  version  of  the  finite  element 
analysis  program  "ANALYZE"  (Ref  18) . The  problem  was  solved  with  a 
mesh  that  consisted  of  452  nodes  and  441  elements.  The  elements  used 
were  quadrilateral  and  triangular  membranes.  Mesh  organization  and 
boimdary  conditions  are  shown  in  Figure  23.  As  a check  of  the  solution, 
displacements  of  the  "pulled"  edge  were  examined.  The  displacements  at 
the  edges  looked  extremely  high  in  comparison  with  displacements  at 
the  center,  but  were  found  to  agree  with  a second  solution  obtained 
through  NASTRAN.  Since  in  our  case  we  have  a very  short  membrane, 
width  = 4,  L ~1.9;  1.27;  1.6;  2.0;  4.0;  and  tQ  = .08,  it  was  suspected 
that  the  restraint  in  the  x direction  has  an  effect  of  singularity. 

To  check  this  point,  the  restraints  in  the  x direction  at  nodes  1, 

441  were  released  and  then  all  restraints  in  x direction,  except  at 
the  mid-node  221,  were  released.  This  modification  did  not  change 
the  results  significantly,  as  can  be  seen  in  Figure  24.  This  suggests 
that  no  singularity  effects  we re  introduced  by  the  x direction  restraint. 
Figure  24  also  shows  the  displacement  variation  as  the  length  of  the 
membrane  grows . 


Appendix  D 


Stepped  Center  Gage  - Stress  Intensity 
Factor  Approximate  Equation 


Figure  25*  Schematic  Stress  Distribution  at  Crack  Tip 

of  Stepped  Gage 


Using  the  relationship  between  stress  intensity  factor  and  stress, 


on  the  line  y = 0f  we  get: 


a 


(69) 


and  tile  force  Pa  produced  over  the  uncracked  region  will  be 


b~a 


o 


Pa 


✓b^a 


*1 
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(71) 


60 


This  force  will  produce  an  average  displacement  on  the  outer  segment 


Uo  - 


P L 
a o 

A E 
o o 


P L 
a o 


The  displacement  of  the  center  segment,  may  be  approximated  by  using 
Rice's  (Ref  22:249)  result  for  an  infinite  strip  with  semi-infinite 


crack, 


Uc  - 


Kx[a  - v2)l c]' 


The  total  displacement  may  be  expressed  as 


Utot  = UC  + uo 


On  the  other  hand,  utot  produced  by  the  displacement  of  the 
bond  line  with  the  carrying  structure  is 


Substituting  and  rearranging  and  making  a nondimensional  stress 
intensity  factor,  normalized  with  respect  to  the  nominal  stress  in  the 


structure  we  get: 
*1 


L + L_ 
o c 
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